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Abstract 

Cross sections for elastic photoproduction of J/ip and T mesons are presented. For 
J j ip mesons the dependence on the photon-proton centre-of-mass energy W yp is analysed 
in an extended range with respect to previous measurements of 26 < W 1P < 285 GeV. 
The measured energy dependence is parameterized as <7 7P oc W^ p with 5 = 0.83 ± 0.07. 
The differential cross section da/dt for J/ip mesons is derived, its dependence on W 7P and 
on t is analysed and the effective trajectory (in terms of Regge theory) is determined to be 
a(t) = (1.27 ± 0.05) + (0.08 ± 0.17) • t/GeV 2 . Models based on perturbative QCD and 
on pomeron exchange are compared to the data. 
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1 Introduction 



A large contribution to the photoproduction of J /if) mesons is the elastic process, jp — > J/if)p, 
which has been studied previously over a wide range, from threshold up to a photon-proton 
centre-of-mass energy, W 1P , of approximately 160 GeV in ep collisions at HERA. The cross 
section for elastic J /if) photoproduction was observed to rise much more steeply with W 1P [[l]-0] 
than that for the light vector mesons [|5p. In perturbative quantum chromodynamics (pQCD), 
the process is viewed as an almost real photon from the lepton coupling to a pair of cc quarks 
which interact with the proton via exchange of two (or more) gluons and evolve into a J /if} 
meson. The cross section is then proportional to the square of the gluon density of the proton. 
The rapid rise of the gluon density at low fractional gluon momenta leads to the observed rise 
of the cross section with W 1P [@-|Sp. In such calculations, the heavy quark mass serves as the 
QCD scale. Therefore, it is of interest to measure the cross section of the T meson, for which 
perturbation theory might be more reliable. In non-perturbative models, J/if> production has 
been described in terms of pomeron exchange [0,[1O||. Since the standard soft pomeron fails 
to describe inclusive electroproduction and the observed W 1V dependence of the J /if) cross 
section, extensions of such models have been proposed JTT]] . In one such extension [|T2|], which 
is compared to the data of this paper, an extra "hard" pomeron is introduced. 

Here, we present new data obtained by the HI experiment which may lead to a better un- 
derstanding of the elastic production process. The W 1P dependence of the cross section for 
elastic J /if) photoproduction is measured, with increased statistics in an enlarged range of 
26 < W JP < 285 GeV, and compared to QCD and pomeron models. Furthermore, the dif- 
ferential cross section da/dt, where \t\ is the squared four-momentum transfer to the proton, is 
measured as a function of t and of W 1P . From its W 1P dependence, an effective Regge trajec- 
tory is determined using data from this experiment alone, thus avoiding relative normalization 
problems when comparing data from different experiments. The results are compared to the 
predictions of the two-pomeron model [ |12| ] and also to a recent calculation in next-to-leading 
order (NLO) using the BFKL equation [|T3|]. Finally, we present a measurement of elastic pho- 
toproduction of T mesons which allows a test of the predictions of pQCD at a higher mass 
scale [|ig,|Bp. 

2 Data Analysis 

We report, here, an analysis of the process ep — > eVp, (V = J /if) or T), where J /if) mesons 
decaying to and e + e~ and /i + pT are observed. The data were taken with the HI 

detector while HERA was operating with positrons of 27.5 GeV and protons of 820 GeV. 

The kinematics of an ep reaction are described by the square of the ep centre-of-mass energy 
s = (p + k) 2 , the squared four-momentum transfer from the positron Q 2 = —q 2 = — (k — k') 2 
and the scaled energy transfer y = (p-q)/(p-k), where k, k', p and q are the four-momenta of the 
incident and scattered positron, the incoming proton and the exchanged photon, respectively. In 
the photoproduction domain Q 2 — > 0, the positron is scattered at small angles and generally 
not observed. In the present analysis, the photoproduction region is defined by Q 2 < 1 GeV 2 . 
The average Q 2 values, derived from the Monte Carlo simulation, are (Q 2 ) ~ 0.05 GeV 2 for 
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the J/ip analyses and (Q 2 ) ~ 0.11 GeV 2 for the T analysis. In the photoproduction limit, the 
square of the centre-of-mass energy of the photon-proton system W 2 p = (p + q) 2 is given by 
W 2 p ~ ys. It is calculated from the relation: y = (E — p z ) v /{2 E e ), where E and p z denote 
the total energy and longitudinal component[] of the momentum of the vector meson V and E e 
is the energy of the incident positron. 

Detector and Data Selection 

The HI detector is described in detail elsewhere [|T^p. The main detector consists of a system 
of drift and proportional chambers: the central track detector (CTD), covering the polar an- 
gle range 15° < 9 < 165° and the forward track detectors (FTD), covering 7° < 9 < 25°. 
The tracking detectors are surrounded by a liquid argon calorimeter (LAr, 4° < 9 < 154°), a 
scintillating fibre calorimeter (spaghetti calorimeter, SpaCal [|T7|], 153° < 9 < 178°) and an 
instrumented iron return yoke (central muon detector, CMD, 4° < 9 < 171°). The interaction 
region is surrounded by an assembly of high precision silicon detectors, the backward part of 
which (backward silicon tracker, BST [|T8"|]) is used in this analysis. The BST, which became 
fully operational in 1997, is a four plane silicon strip detector telescope arranged concentrically 
around the beam axis to measure the polar angle of tracks (172° <9< 176.5°). 

Electrons are identified in the LAr or the SpaCal calorimeters. Muons are identified in the LAr 
calorimeter, in the CMD and, at small polar angles, in the forward muon detector (FMD, 3° < 
9 < 17°), which consists of drift chambers on either side of a toroidal magnet. The momentum 
and angular measurements of the decay leptons are provided by the CTD (FTD) in the central 
(forward) detector region. At the highest values of W 1P , both electrons are scattered into SpaCal. 
In this case, the polar angle and the event vertex are derived from track measurements in the 
BST. 

In elastic photoproduction of J/ip and T mesons, both the scattered positron and the scattered 
proton are not generally detected. The largest background is due to processes in which the 
proton dissociates to form a hadronic state Y, the proton remnant. To reject such events, the so- 
called "forward detectors" are utilized: the proton remnant tagger, PRT, scintillation counters 
covering polar angles close to the proton beam, the drift chambers of the FMD in front of the 
toroid and the forward part of the LAr calorimeter (9 < 10°). Proton remnants with masses 
My > 1.6 GeV can be detected. 

The experimental signature of the process ep — > eV p, V — > £ + £~ in the photoproduction do- 
main is a pair of leptons in an otherwise "empty" detector. The data selection thus requires only 
two decay leptons to be observed in the main detector with no signal in the forward detectors. 
The value of W 1P is related to the polar angles of the vector meson and therefore also to the po- 
lar angle of the decay leptons. Selection criteria are used which correspond to different detector 
regions and ranges of W 1P (see Table [I]). For each analysis region, the decay mode, or 
e + e~, with the better detection efficiency is chosen. 

In the central region, J/ip — > is analysed in the range 40 < W yp < 150 GeV. Exactly 

two oppositely charged particles, measured in the CTD, are required with transverse momenta 

'The forward (+z) direction, with respect to which polar angles are measured, is defined as that of the incident 
proton beam direction. 
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Table 1: Summary of the different data sets. 



(with respect to the beamline) p t > 0.8 GeV. At least one of these must be validated as a 
muon in the LAr calorimeter or in the CMD (for further details see [|l],|3|, [19|]). Background 
from cosmic ray muons is rejected by an acolinearity cut. Triggers based on muon and track 
signatures from the decay leptons are used. For the reconstruction of T decays, the acceptance 
region is 70 < W 1P < 250 GeV. In order to optimize the signal to background ratio, both 
muons must be identified and a tighter cut against cosmic rays is applied. 

For the J/ ip analysis in the backward region, the decay to electrons is used. Two sets of selection 
criteria are applied, depending on the event kinematics. The first set (135 < W 1P < 210 GeV) 
requires one decay electron to be measured in the CTD and one in the SpaCal calorimeter. 
The electron measured in the CTD must have a momentum p > 0.8 GeV and a transverse 
momentum p t > 0.7 GeV. It is validated by an electromagnetic cluster in the LAr calorimeter. 
The other electron is selected by requiring a cluster in the SpaCal calorimeter with energy 
above 4.2 GeV, which is assumed to originate from the intersection of the first electron track 
with the z-axis. In order to suppress inelastic contributions and background, no further tracks 
are allowed and any energy in the SpaCal outside the selected electron cluster must be below 
6 % of its energy. The second selection (210 < W 1P < 285 GeV) requires both electrons to be 
detected in the SpaCal calorimeter as described above. They must be validated either by hits in 
a proportional chamber or by a track in the BST By requiring two charged particles, at least one 
of which is in the acceptance of the BST, most of the non-resonant background from the QED 
process ep — > ejp is rejected. The triggers for both selections are based on signals from the 
SpaCal and the CTD. In addition the triggers in the central and backward regions use a trigger 
based on neural networks [ |20| ] . 

In the J/ip analysis in the forward region (26 < W rp < 36 GeV), the FMD is used to identify 
one decay muon with momentum above 5 GeV. The other muon is measured in the CTD 
with momentum above 0.8 GeV. No other tracks, except those associated with the muons, 
are allowed. Proton dissociative events are rejected by a forward detector analysis, taking into 
account that one of the muons passes through the drift chambers of the FMD and may pass 
through the forward part of the LAr calorimeter. The FMD also supplies the trigger for these 
events in conjunction with signals from track chambers and the central muon detector. 
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Figure 1: Mass spectra for the four W 1P regions of the elastic J/tp selection: a) pairs in 
the central region, b) and c) e + e~ pairs in the backward region, d) pairs in the forward 
region. A fit of the signal region as described in the text (full line) and the simulated mass 
spectra from QED background processes (shaded histograms) are shown. 

Cross Section Determination 



After the data selection, clear peaks at the J/tp mass are observed in the distributions of the 
invariant mass of the two leptons in all analysis regions (Figure p. The remaining background 
in the forward and central regions is dominated by the process 77 — > where the photons 

originate from both the positron and proton. In the backward region at very high W 1P , the 
background is composed of the processes 77 — > e + e~ and ep — > cyp, where the photon and 
electron are scattered at large angles. These background processes can be simulated by the 
LPAIR [JHJ] and COMPTON [|Z2p generators, respectively. 

The number of J/tp events in the central region (Figure [l]a) is determined in each analysis bin 
by a fit which uses a Gaussian for the signal and a polynomial for the background. In the 
backward region (J/tp— > e + e~), the number of events is determined by fitting the signal region 
with a Gaussian and an exponential to describe the low mass tail, and subtracting the simulated 
background (Figure |l]b,c). The low mass tail is caused by radiative effects in the the detector 
material. In the forward region, the background is determined by the Monte Carlo simulation. 

For the calculation of the cross section, the number of events without a signal in the forward 
detectors is corrected for acceptance and efficiency losses by using the Monte Carlo simulation 
DIFFVM [E3p. DIFFVM is based on the Vector Meson Dominance Model and Regge phe- 
nomenology. Either the elastic or the proton dissociative processes can be generated. The Q 2 , 
W 1P and t dependences have been adjusted to fit the data. 

A correction is applied for the remaining background from proton dissociation which is typi- 
cally 10 — 13% in the central and backward regions and 30% in the forward region. It is deter- 
mined from the Monte Carlo simulation and cross checked with data samples with and without 
signals in the forward detectors. The non-prompt J/tp events from tp{2S) decays (~ 3%) are 
subtracted according to their measured cross section [ |24| ] and branching ratio [[25]]. The total 
efficiencies (including acceptance) for the J/tp analyses are typically 5%, 14% and 20% in the 
forward, central and backward regions, respectively. The trigger efficiencies are about 18%, 
20% and 38 - 60%, depending on W 1P . At the highest W 1P values (W JP > 200 GeV), the 
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trigger efficiency is 100%, within errors. The losses in this region are mainly due to the limited 
acceptance. 

The ep cross section is obtained using the measured luminosity and the branching ratio of the 
J/ip decay to electrons or muons, (6.0 ± 0.2)% [ |2~5| ] each. For the decay to electrons, the 
radiative decay J/ip — > e + e~7 with the branching ratio (0.88 ± 0.14)% [|25|] is also taken 
into account. The photoproduction cross section is calculated from the electroproduction cross 
section, assuming factorization of the ep reaction into emission of photons, described by a 
photon flux (Weizsacker- Williams approximation P5[]), followed by a ^p reaction. 

The systematic error for the J/ip analysis in the central region increases with W 1P from 11% 
to 14% with an estimated correlated part, which affects only the normalization, of 5 — 6%. 
The largest contributions are the uncertainties in the determination of the proton dissociative 
background (5—10%), the efficiency of the lepton identification (6.5%) and the trigger efficiency 
(6%). In the backward region, the systematic uncertainties are estimated to be, on average, 16%. 
For W 1P < 200 GeV, the contributions are similar to those in the central region. For W 1P > 
200 GeV, the main contributions are due to the uncertainties in the subtraction of the residual 
background from QED processes (10%) and proton dissociation (10%). In the forward region 
(low W 1P ), the total uncertainty is 26%, where the largest contributions are the uncertainty of 
the trigger efficiency (18%) and the determination of the background (16%). 

Radiative corrections have been neglected throughout the analysis. Recent calculations for 
elastic photoproduction of p mesons [ |2"7| ] suggest that, for J/ip photoproduction, they are small 
compared to experimental errors. 



3 W 1P Dependence of J/jJj Production 

The cross section for elastic photoproduction of J/ip mesons, as a function of W 1P , is listed in 
Table ^ and is shown in Figure ||a, together with previous HI and ZEUS results [|I],|2|]. The data 
show good agreement in the region of overlap. A fit of the form W^ p to the present HI data 
between 26 and 285 GeV yields a value of 5 = 0.83 ± 0.07 (including statistical and systematic 
errors). This result confirms, with smaller errors, the steep rise of the J/ip cross section with 
W 1P observed previously [[]],§], which is much larger than for the light vector mesons [§]. 

In Figure ^a, the HERA data are shown with predictions from a leading order pQCD calcula- 
tion [0,0] using various gluon density distributions. The important prediction concerns the W 1P 
dependence of the cross section since the absolute magnitude is more dependent on a number of 
parameters and on non-perturbative effects. The slope of the data is described well using either 
the gluon density CTEQ4M ^ or MRSR2 @ while the GRV-HO ^ based result is too 
steep. 

For the comparison with models based on pomeron exchange, the present data, previous HERA 
data and results from fixed-target experiments [|3T]]Q are used (Figure ^b). A pomeron trajec- 
tory of the form a(t) = a + a't leads to a W 1P dependence o 1P oc W^ p ao ~^ /b(a' \ \nW lp ) 



2 Only fixed-target results are shown which were performed on H2 or D2 targets and which have been corrected 
for contributions from proton dissociative processes. 
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W TD [GeV] W ra [GeV] 



Figure 2: a) 77ze cross section <j(-yp — > J/ V>p) versus W 1P from this analysis and previous HERA 
results. The inner and outer error bars show the statistical and total errors, respectively. The 
full line represents a fit of the form <r 7P oc W* to the present HI data, yielding 5 = 0.83 ± 0.07 
with a x 2 / dof = 4.6/13. QCD predictions by Frankfurt et al. /f|/ using various parameteriza- 
tions of the gluon density in the proton and the charm quark mass of 1. 48 GeV as chosen in /^/ 
are also shown, b) The same HERA data as in a) and results from fixed-target experiments with 
a fit of the two-pomeron model by Donnachie and Lands hoff {full line; \ 2 ' j 'dof = 16.1/22). 
Only the contributions of the soft and hard pomeron amplitudes were adjusted allowing for 
mixing /|72|/. The separate contributions are indicated. 

[J3|.[r0p. The parameters of the established soft and the proposed hard pomeron trajectories are 

Oo, a') S oft = ( L08 > °- 25 GeV ~ 2 ) and («o, «')hard = (1-418, 0.1 GeV~ 2 ) [0. By using these 
trajectories and fitting the relative contributions and mixing as defined in JT^j], a good represen- 
tation of the data is given. The relative contributions of hard, soft and mixing terms are found 
to vary between 0.1 : 0.5 : 0.4 at W yp = 30 GeV and 0.5 : 0.1 : 0.4 at W yp = 250 GeV. 

4 Differential Cross Section da / dt for J/i/j Mesons 

The dependence of the photoproduction cross section of J/tp mesons on t is studied in the 
region of 40 < W 1P < 150 GeV, which provides good t resolution. First, this dependence is 
measured averaging over the entire W 1P range, assuming an exponential behaviour. Then, the 
dependence of da/dt on W 1P in fixed t bins is determined and an effective Regge trajectory is 
derived. 

t Distribution: The variable t is approximated by t ~ —(p t ) 2 , where p t is the transverse 
momentum of the J / ip meson. The resolution of t, obtained from the Monte Carlo simulation 
of vector mesons, is typically 0.035 GeV 2 and increases to 0.060 GeV 2 for \t\ w 1.2 GeV 2 . 
The events without a signal in the forward detectors, in the mass range 2.9 < M w < 3.3 GeV, 



9 



£io 2 



10 



i — 1 — 1 — 1 — i — 1 — 1 — 1 — i — 1 — 1 — 1 — r 


— 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — r 

• H1 




-- fit 




— elastic 




pdiss. 




non-res. 


! . . , 1 , , , 1 '\. , , 


\ v ~ • 



0.2 



0.4 



0.6 



0.8 



1.2 

2i 



lit I [GeV z J 



Figure 3: The differential cross section da j dtfor elastic J/ ip photoproduction, as a function of 
\t\, averaged over the energy region 40 < W 1P < 150 GeV. Only statistical errors are given. 
Also shown is the result of a fit (dashed line; y 2 / dof = 4.7/6) to the data which includes the 
three different contributions described in the text. The solid curve shows the fit result for elastic 
J j ijj production assuming an exponential t dependence. 



are divided into bins in t and corrected for efficiencies but not for remaining backgrounds. 
Photon-proton cross sections are calculated as described in section |TJ, yielding the resulting 
values of da/dt shown in Figure ^|. 

With this selection, the data contain elastic J/ip events (81%), events with proton dissociation 
(13%) and non-resonant background (6%). These fractions result from the analysis described in 
section |T] and are fixed in the subsequent analysis of the t-dependence. The t-dependence of the 
proton dissociative J/ip contribution is obtained by fitting those corrected data which contain 
a signal in the forward detectors with one exponential function. The different t dependence of 
the efficiencies for proton dissociative events with and without signals in the forward detectors 
is taken into account. The non-resonant background is described by a sum of two exponential 
functions which are obtained by fitting the data outside of the J/ip mass region. 

A common fit of the three contributions to the data is then carried out, assuming an exponential 
distribution e~ b ''! for the elastic cross section (see Figure 0). The only free parameters in this 
fit are the elastic slope parameter b and an overall normalization. Only the statistical errors of 
the data are taken into account in the fit. A value b = (4.73 ± 0.25+Q3g) GeV~ 2 is derived for 
\t\ < 1.2 GeV 2 , where the first error is statistical. The systematic uncertainty is obtained by 
varying the data selection cuts, the t distribution of the generated events for both the acceptance 
and efficiency corrections and the admixture of the background contributions (both relative 
normalizations and shapes). The resulting value of b agrees within errors with the previous HI 
result [[]]] and with a similar result from the ZEUS collaboration [0. 
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Regge Trajectory: In the description of elastic scattering based on Regge phenomenology 
and pomeron exchange, the energy dependence of the elastic cross section follows a power law: 



da 



da 
~dt 



-bo |*| 



w. 



IP 



w 



4(a(t)-l) 



(1) 



where a(t) = ao + a't denotes the exchanged trajectory and b and Wo are constant parameters 
not determined by the theory. Whereas Regge theory predicts a change of the W 1P dependence 
with t, there is no such simple relationship in pQCD 0. 

In order to determine the effective trajectory a(t) in elastic J/ip production, the data used for 
Figure ^| are analysed as functions of W 1P and t. The data are divided into five bins in t and 
six bins in W 1P . According to the Monte Carlo simulation, the migrations due to resolution 
effects range from 20 — 50% in t and are negligible in W JP . The differential cross section da/dt 
is determined in each bin (tj, Wf), as described in section |lj, i.e. correcting for non resonant 
and proton dissociative background. The resulting values are listed in Table [3] and shown in 
Figure ^ as functions of W 1P in the five t bins. Statistical and systematic errors, excluding the 
contributions which affect only the normalization, are taken into account; the statistical errors 
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Figure 4: The differential cross section da/dt as a function ofW lp in five bins oft together with 
a fit of the form da/dt = N ■ (W rp /Wo)^ a ^^ (solid line). The inner error bars on the data 
points show the dominating statistical errors and the outer bars the total error. The predictions 
of the soft and hard Donnachie-Landshoff pomeron trajectories /[72[/ are shown as dash-dotted 
and dashed lines. The prediction based on a NLO BFKL calculation /|73|/ is given as a dotted 
line. All theoretical curves are normalized to the fit at Wo = 90 GeV. 
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Figure 5: The measured Regge trajectory a(t) = «o + a't for the process jp — > J/ipp. The 
solid line shows the result of the fit. The one standard deviation contour is indicated by a shaded 
band. Also shown are the soft and the hard Donnachie-Landshoff pomeron trajectories /|72[/ and 
a result based on a NLO BFKL calculation /[73[/. 

dominate. The data in each bin U are fitted as a function of W 1P (Equation [j]) with a(U) and a 
normalization constant as free parameters. 

The results of these fits are shown as full lines in Figure^. The respective curves, corresponding 
to the soft and the hard pomeron trajectories by Donnachie and Landshoff [ |12| ] as well as a 
hard pomeron extracted from a NLO BFKL calculation [O], are also shown. Note that the 
differential cross section da/dt at fixed t exhibits a rise with W 1P at all t values, in contrast to 
the expectation from the soft pomeron model. 

The resulting fit values of a(t«) are shown as points in Figure [5|; the error bars contain the 
statistical and systematic uncertainties. A fit to these values of the form a(t) = a + a't, where 
a and a' are free parameters, yields the result 

a(t) = (1.27 ± 0.05) + (0.08 ± 0.17) • t/GeV 2 . 

This fit is shown in Figure [5|, including a band which reflects the one standard deviation un- 
certainty, taking into account correlations between a an d a'. The resulting effective trajectory 
lies, as expected from the analysis of the total cross section (see Figure [Z]a), in between the soft 
and the hard pomeron trajectories of Donnachie and Landshoff. A prediction for a pomeron, 
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derived from a NLO BFKL calculation [[T3|], is also shown in Figure H. Note, however, that in 
this calculation only a is predicted^ and a' = is assumed for the plot. The measured inter- 
cept a in J/ tp production is found to be significantly larger than in the production of the light 
vector mesons p and </>. The slope a' is found to be compatible with zero, in contrast to the value 
measured for p and </> mesons [^2p. The present results allow, for the first time, the extraction 
of the effective Regge trajectory exchanged in J/ip photoproduction from a single experiment, 
avoiding relative normalization uncertainties between experiments. A previous fit [B3J], using 
earlier HERA results and fixed target data at lower W 1P , gives compatible results. 



5 T Production 

The mass distribution of the selected di-muon events, without a signal in the forward coun- 
ters, is shown in Figure |6|a. An excess of events is visible at a mass of ~ 9.4 GeV above a 
background which is fitted by a power law. The background is, within errors, described by 
the process 77 — > p + p~. The significance of the excess is about two standard deviations 
and its position is compatible with these events arising from T(15) decay. The mass resolu- 
tion is a ~ 250 MeV. Due to limited statistics, small contributions from the two lowest radial 
excitations T(25) and T(35) cannot be excluded. In order to estimate the number of signal 
events, a wide mass interval is chosen which extends from 2 a below the T(IS') to 2 a above 
the T(3S) state (8.9 < m m < 10.8 GeV). Subtracting the background, 12.2 ± 6.3 signal events 
are obtained. Corrections for the remaining background, due to diffractive proton dissociation 
(16%) and the total selection and trigger efficiency (24%), are determined from the Monte Carlo 
simulation, cross checked with data similarly to the J/ip analysis. Taking into account the total 



integrated luminosity and using the Weizsacker- Williams approximation [26], a 7^ cross section 



in the range 70 < W 1P < 250 GeV with an average value of (W 1P ) = 143 GeV is obtained: 
<r(7 p -> T p) x BR(T -> p + p~) = (19.2 ± 9.9 ± 4.8) pb, 



where T includes the states T(IS'), T(2S) and T(3S). The first error is statistical and the second 
systematic (25%). The latter is calculated by varying the definition of the signal and background 
regions, using different estimates of the background function or calculating the background 
from the 77 — > p + p~ prediction and taking into account uncertainties of efficiencies. 

In order to extract the cross section for the production of the T(IS') state, an assumption must 
be made about its relative production ratio. We choose a fraction of 70%. This is in broad 
agreement with an estimate using the branching ratios and the electronic widths of the different 
T states [ |25| ] and with a recent calculation [pfl. Using the standard branching ratios for T(IS'), 



T(2S) and T(3S') Q25|], the result for the elastic T(IS') photoproduction cross section is 

a(7p -> T(1S» = (0.55 ± 0.28 ± 0.14) nb . 

In Figure |^b), the measured cross section a{^p — > T{lS)p) is shown together with the ZEUS 
measurement [[34]], which agrees well with our result. Recent pQCD calculations [|T3],[r5|] are 

3 We chose the solution corresponding to £ = 3 in Table 2 of Jl3|]. 
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able to describe the data after consideration of two effects: the real part of the scattering ampli- 
tude and the non-diagonal parton distributions in the proton. Such effects are found to be more 
important for the production of the T than for the J /if) meson, due to its larger mass. 



6 Summary 

Improved results over an extended kinematic range on the elastic photoproduction of J/ip 
mesons are presented. The total photoproduction cross section is measured as a function of 
the photon-proton centre-of-mass energy W 1P in the range 26 < W 1P < 285 GeV, where it can 
be parameterized Wfp with 5 = 0.83 ± 0.07. 

The data are well described by a perturbative QCD calculation by Frankfurt et al. [|8|], where a 
good description is obtained using the CTEQ4M or MRSR2 gluon densities in the proton. The 
data in this W 1P range can also be described by a non-perturbative model which assumes the 
exchange of two pomeron trajectories. However, a model with only a soft pomeron is ruled out. 



The differential cross section, da/dt, for elastic J ftp production, averaged over the range 40 < 
W JP < 150 GeV, is well described by a single exponential for the accessible range below 
\t\ = 1.2 GeV 2 with a slope parameter b = (4.73 ± 0.25±g;^) GeV~ 2 . 




[GeV] W rp [GeV] 
Figure 6: a) The invariant mass spectrum of the selected di-muon events without a signal in 
the forward detectors. The solid line shows the result of a fit to the background region. The 
histograms are Monte Carlo predictions for 77 — > fi + pr(LPAIR) and for the three lowest T 
states (DIFFVM, normalized to the observed number of events). The arrows indicate the chosen 
signal region, b) The cross section for elastic photoproduction of T (IS) from this analysis with 
the ZEUS result (which also assumes a 70% contribution of the IS state). The inner and outer 
error bars show the statistical and the total errors, respectively. Predictions based on pQCD, 
calculated with the gluon distribution MRSR2 /j2"P|/, are shown. MRT( 1 ) /]75|/ and FMS /|7^[/ 
are calculations based on the leading order vector meson cross section, including corrections. 
MRT(2) /|75[/ employs parton hadron duality to derive the prediction for the T from the bb cross 
sections. 
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In the same W 1P range, da/dt is measured as a function of W JP and t. From this analysis, the 
parameters of the effective trajectory for elastic J/ip photoproduction are determined using data 
from this experiment only. The values of the intercept and the slope of the measured trajectory 
lie between those of the well known soft and the conjectured hard pomeron trajectories of 
Donnachie and Landshoff. The measured intercept is compatible with a prediction derived 
from a NLO BFKL calculation. 

Elastic photoproduction of T mesons is observed. The cross section for the sum of the three 
lowest states is measured. The elastic photoproduction cross section for T(IS') mesons is ex- 
tracted. Recent calculations within pQCD are in agreement with this value. 
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interval [GeV] 


(W w ) [GeV] 


(7(7^ 


- J/#) [nb] 


26 


- 36 


31.0 


19.8 


±3.7 


±5.1 


40 


- 52 


46.2 


29.8 


±3.7 


±3.4 


52 


- 60 


56.1 


41.7 


±5.2 


±4.8 


60 


- 67.8 


64.0 


49.5 


±6.6 


±5.6 


67.8 


- 77 


72.5 


51.8 


±6.9 


±5.9 


77 


- 86 


81.5 


62.4 


±8.5 


±7.1 


86 


- 97 


92.0 


67.6 


±8.8 


±7.7 


97 


- 113 


105.2 


64.6 


±9.6 


±7.4 


113 


- 150 


133.4 


89.0 


± 13.2 


± 12.7 


135 


- 160 


147.3 


80.6 


±7.9 


± 12.7 


160 


- 185 


172.2 


104.0 


±7.9 


± 16.4 


185 


- 210 


197.1 


113.6 


± 12.2 


± 18.0 


210 


- 235 


222.4 


116.0 


± 14.8 


± 18.5 


235 


- 260 


247.4 


118.9 


± 15.6 


± 19.0 


260 


- 285 


272.4 


156.9 


±20.3 


±25.0 



Table 2: Cross sections for the elastic process 7p — > J / 'tp p in bins ofW lp . The first error of the 
cross section is the statistical error and the second is the systematic uncertainty. 









W rp interval [GeV] 




-t [GeV 2 ] 


40-54 


54-66 


66-78 


78-91 


91 - 110 


110-150 


0-0.053 


146.7 


231.6 


392.4 


220.1 


348.2 


553.5 


(0.03) 


±32.5 


±48.4 


±89.3 


±58.7 


±90.1 


±138.7 




±13.4 


±21.1 


±35.7 


±20.0 


±31.7 


±58.1 


0.053 - 0.134 


103.8 


102.6 


140.7 


224.9 


202.6 


219.4 


(0.09) 


±21.3 


±20.3 


±80.9 


±45.2 


±43.3 


±64.9 




±9.4 


±9.3 


±12.8 


±20.5 


±18.4 


±23.0 


0.134-0.25 


48.9 


62.0 


80.9 


125.2 


93.4 


147.6 


(0.19) 


±10.9 


±13.6 


±18.6 


±25.0 


±23.3 


±38.9 




±4.4 


±5.6 


±7.4 


±11.4 


±8.5 


±15.5 


0.25 - 0.44 


41.3 


40.0 


29.5 


57.2 


64.6 


94.6 


(0.34) 


±8.8 


±9.4 


±7.5 


±12.6 


±16.0 


±23.1 




±3.8 


±3.6 


±2.7 


±5.2 


±5.9 


±9.9 


0.44- 1.2 


4.82 


3.82 


5.42 


9.62 


5.72 


11.43 


(0.72) 


±1.09 


±0.99 


±1.33 


±2.54 


±1.81 


±3.35 




±0.44 


±0.35 


±0.49 


±0.88 


±0.52 


±1.20 



Table 3 : Differential cross sections da/dt[ nb/GeV 2 ] for the elastic process •yp — > J/ipp in bins 
of t and W 1P . The t intervals and bin centres are given in the left column. The first error is 
statistical and the second is the systematic uncertainty neglecting pure normalization errors. 
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